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Introduction 

The  reduced  cohesion  of  grain  boundaries  (GBs)  is  often  the  controlling  factor  limiting  ductility, 
and  hence  performance  and  reliability  of  high-strength  metallic  alloys  (1).  Intergranular 
embrittlement  in  metals  is  usually  caused  by  impurities  segregating  towards  the  GBs  (2-6).  A 
dtictile-biittle  transition  temperature  (DBTT)  as  low  as  -196C  (7)  was  observed  in  high  purity  W  single 
crystals  obtained  by  electron  beam  zone  melting  with  special  impurity  gettering.  Impurities  present 
in  bulk  concentrations  of  lO'^-lO*^  atomic  percent  can  result  in  a  dramatic  decrease  of  plasticity, 
drastically  degrading  mechanical  properties  of  metallic  alloys,  and  thus  posing  significant 
tecbmological  and  application  problems.  This  detrimental  effect  of  minute  impurity  concentrations 
can  be  readily  understood.  A  simple  estimate  shows  that  a  ppm  amount  of  impurity  is  sufficient  for 
saturating  aU  the  grain  boundaries  in  a  typical  grain-size  polycrystal.  Sensitivity  of  the  DBTT  to  the 
grain  size  confirms  the  above  physical  concept:  the  larger  the  grain  size,  the  smaller  amount  of 
impurity  is  needed  to  saturate  the  GB  (3).  Fine-grain  polycrystals  are  known  to  be  less  brittle.  It 
should  be  noted  that  BCC  crystals,  being  not  as  dose  packed  as  FCC  or  HOP,  are  partictilarly  prone 
to  GB  embrittlement  by  impurity  segregation. 

If  impurities  are  the  main  cause  of  embrittlement,  gettering  the  "the  harmful"  impurities, 
such  as  O,  N,  P,  S,  As,  Sb,  etc.,  from  the  GBs  by  their  forming  thermodynamically  stable  phases  with 
other  elements,  is  an  obvious  way  of  improving  ductility.  This  process,  however,  requires  careful 
control  since  ^e  ductility  upon  gettering  will  be  improved  oi^y  so  far  as  these  second  phase 
predpitates  remain  fine;  any  coa^ation  of  predpitates,  such  as  the  so-called  Ostwald  ripening, 
would  result  in  an  adverse  embrittling  effect . 

Fe-base  alloys  and  W  are  among  the  most  important  ultra-high-strength  materials.  The 
purpose  of  this  paper  is  to  eluddate  ffie  energetics  of  impurities  in  Fe  and  W  GBs,  and  analyze  the 
effect  of  impiurities  on  the  intergranular  cohesion  at  the  electron-atom  level.  A  deeper  understanding 
of  the  cohesion-decohesion  processes  on  the  microscopic  level  will  lay  a  foundation  for  a  "smart 
design"  of  more  ductile  and  tough  alloys.  In  particular,  our  theoretical  analysis  of  the  electron 
structure  and  the  energetics  of  Fe  and  W  GBs,  both  clean  (CL)  and  with  impurities,  enables  one  to 
make  important  predictions.  Our  calculations  have  shown  that  impurities,  such  as  N,  O,  P,  S,  Si 
weaken  the  intergranular  cohesion  resulting  in  "loosening"  of  the  GBs  in  Fe  and  W.  The  presence  of 
B  on  the  contrary,  enhances  the  interatomic  interaction  across  the  GB.  The  so-called  "site- 
competition  effect"  should  play  an  important  role  affecting  impurity  distribution  in  GBs.  Among  the 
impurities  analyzed,  both  in  Fe  and  W,  B  in  the  GB  has  the  lowest  energy,  and  thus  would  tend  to 
push  other  impurity  atoms  off  the  GB,  ffius  improving  resistance  to  brittle  fracture. 

Energetics  of  Imnurities  in  Fe  and  W  GBs  and  the  Site  Competition  Effect 

The  progress  in  developing  efficient  methods  of  first-principles  calculations  and 
computational  algorithms  have  made  possible  systematic  studies  of  the  role  of  impurities  in 
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intergranulai*  cohesion  of  transition  metals  at  the  atomic  and  the  electron-ion  level.  Calculations  on 
both  cluster,  two-dimensional  and  supercell  models  of  GBs  with  impurities  have  provided  an  in- 
depth  insight  into  mechanisms  of  GB  decohesion  (for  references,  see  (8,9)). 

Since  the  first-principles  electronic  calculations  on  low-symmetry  systems  (such  as  lattice 
defects  or  GBs)  are  still  extremely  complicated  and  costly,  semi-empirical  methods  based  on  reliable 
first-principles  foundations  have  also  been  developed.  Among  them,  the  most  popular  is  the 
Embedded  Atom  Method  (EAM)  (10,11).  This  method  has  been  successfully  used  in  a  wide  variety  of 
calculations. 

In  order  to  study  the  energetics  of  impurity  atoms 
in  Fe  and  W  GBs,  we  have  chosen  first  to  calculate  the 
quantity  which  may  be  called  the  "environment-sensitive 
embedchng  energy"  (ESE)  -  the  energy  of  an  impurity 
atom  in  an  atomic  environment  typical  for  a  GB. 
Knowledge  of  these  energies  for  various  impurities 
enables  one  to  compare  the  relative  stability  of  a 
particular  impurity  in  the  GB  environment. 

Having  calculated  the  ESEs  for  a  number  of 
impurity  atoms,  one  can  use  this  information  in  a 
modified  EAM  approach  for  calculating  the  GB 
relaxation.  The  latter  calculation  enables  one  to  draw 
important  conclusions  regarding  the  intergranular 
cohesion  in  Fe  and  W  in  the  presence  of  a  specific 
impurity  in  the  GB. 

The  model  chosen  for  the  GB  environment  is  an  8- 
atom  hexagonal  supercell  (Me0X,  where  Me  is  Fe  or  W, 
and  X  is  an  impurity  atom).  The  supercell  is  shown  in 
Fig.l,  together  with  the  capped  trigonal  prism  coordi¬ 
nation  of  the  surrounding  Me  atoms. 

A  tri^nal  prism  GB  configuration  is  believed  to  be  a  typical  GB  environment  in  BCC  metals 
and  is  pre<ficted  by  the  theory  of  hard  sphere  packing.  Atomistic  relaxation  studies  have  shown  that 
in  Fe  an  impurity  atom,  such  as  P  or  B,  is  likely  to  occupy  an  interstitial  position  in  the  center  of  ^e 
trigonal  prism  formed  by  Fe  atoms  in  the  GB  core  (even  if,  like  in  the  case  of  P  and  B,  the  impurity 
forms  a  substitutional  solid  solution  with  the  host).  The  hexagonal  supercell  of  Fig.  1  has  a  relatively 
hi^  syi^etry;  it  also  emulates  a  (111)  £3  GB  environment  (8).  This  model,  however,  does  not  allow 
segregation  of  tii^Uy  bound  pairs  of  different  impurity  atoms. 

We  performed  the  spin-polarized  scalar-relativistic  Linear  Muffin  Tin  Orbitals  (LMTO) 
calculations  (our  method  and  approximations  were  the  same  as  in  (8)).  First,  a  series  of  calculations 
(for  six  different  volumes)  were  performed  with  an  impurity  absent  from  the  supercell,  i.e.  an  empty 
sphere  of  the  same  radius  as  that  of  the  radius  of  the  impurity’s  atomic  sphere  was  substituted  for 
the  latter.  Similar  calculations  were  then  performed  for  each  of  the  impurities:  H  (Fe  GB),  B,  C,  N, 
O,  Al,  Si,  P  and  S.  The  ESEs  were  defined  as  follows; 

ESE  =  E(Me6«)  •E(Me60)-E(«)  (1) 

where  E(Meg#)  and  E(MegO)  are  respectively  the  energies  of  the  supercell  with  and  without  the 
impurity  (O  stands  for  an  empty  sphere  substituted  for  the  impiuity  atom),  and  E(#)  is  the  energy  of 
the  firee  impurity  atom.  In  order  to  make  the  calculations  more  consistent,  we  have  chosen  to  use,  as 
E(9)8,  the  values  of  E(Meg€)  -  ElMegO)  extrapolated  to  the  zero  charge  density  (nsO),  which  would 
correspond  to  the  energies  of  impurities  in  the  GB  environment  with  the  host  crystal  lattice  infinitely 
expanded.  The  ESE  energies,  Eq.(l),  as  a  function  of  n,  the  electron  charge  density  due  to  the  host 
atoms  at  the  impurity  site,  are  presented  in  Fig.  2. 

Plots  in  Fig.  2  explain  an  experimentally  observed  phenomenon  known  as  the  "site  competition 
effect ".  As  one  can  see,  in  the  range  of  electron  charge  density  typical  of  a  Fe  GB  (0.015-0.025  a.u.),  B 
and  C  have  the  lowest  ene^  and  thus  would  tend  to  displace  ^e  other  impurities  off  the  GB.  Thus, 
there  exists  a  "site  competition  hierarchy".  In  fact,  in  Fe,  C  is  known  to  successfully  compete  with  P 


Hg.  l.The  Me^X  hexagonal  supeicell  eniulating  a 
typical  trigonal  prism  environment  of  Fe  or  W 
aioms.in  the  (1 1 1)£3  GB:  a)  the  supeicell;  b)  the 
trigonal  prism  cocxdination;  ®  Me,  •  Impurity. 
Numben  designate  atoms  Mel  through  Me4 
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(12)  and  S  (13,14),  and  N  --  with  P  (IS)  and  S  (13).  B  should  obviously  be  more  efficient  in  pushing  the 
harmful  impurities  off  the  GB.  It  was  found  recently  (16,17)  that  B  successfully  competes  with  P. 
Since,  in  principle,  ESEs  are  to  be  sensitive  to  the  type  and  orientation  of  the  GB,  in  some  GBs  B  and 
C  may  compete  in  a  reverse  way.  Also,  it  can  be  seen  that  H  would  successfully  compete  with  all  the 
impurities  but  B  and  C,  N  being  a  maiginal  case. 

In  W  B  has  the  lowest  energy  among  the  impurities  analyzed,  and  the  "gap"  between  B  and  C 
is  significant.  There  is  also  a  site-competition  hierarchy.  It  was  confirmed  experimentally:  N  was 
found  to  successfully  compete  with  C  (18);  while  C  competes  with  P  (19). 


Fig.  2.  The  "environment  sensitive  embedding  energies,  ESE,  vs  n,  the  electron  charge  density 

(in  atomic  units) 


The  plots  in  Fig.  2  also  reveal  an  important  aspect  of  GB  impurity  behavior.  All  the  plots  have 
well  pronounced  minima.  The  positions  of  the  minima  correspond  to  the  electron  density  at  the 
impurity  site  due  to  the  surrounding  Me  atoms  which  would  occur  if  the  GB  were  allowed  to  relax  in 
such  a  way  as  to  minimize  the  impurity's  energy.  The  minima  positions  systematically  (except  for  H 
in  Fe)  shift  towards  lower  densities  with  the  impurity  losing  its  competitive  power.  A  smaller  charge 
density  means  a  "looser"  GB,  less  strong  and  more  prone  to  decohesion.  The  minimization  of  the 
total  GB  energy  (rather  than  only  the  energy  of  the  impurity  atom)  gives  the  characteristic  charge 
densities  which  are  somewhat  higher  than  those  in  the  minima.  However,  from  this  point  of  view, 
N,  0,  S,  P,  Si  and  Al  are  the  obvious  candidates  for  being  "decohesive",  N  is  marginal,  while  B  and  C 
may  be  called  "cohesion  enhancers".  In  fact,  B,  C  and  N  were  experimentally  found  to  improve  the 
GB  cohesion  in  both  Fe  (12-17)  and  W  (18-22),  while  O,  Si,  P  and  S,  being  strong  embrittlers  both  in  Fe 
and  W  (3,6.23,24),  are  believed  to  weaken  the  GB  cohesion  . 

As  mentioned  above,  the  GB  environment  we  were  dealing  with  was  that  of  the  (111)  I3  tilt  GB. 
The  GB  structure  can  be  represented  as  a  succession  of  (111)  hexagonal  planes; 


....CBACBACBACBACBABCABCABCABCABC.... 


(the  GB  plane  is  marked  by  A).  The  CBABC  atomic  stacking  of  the  core  of  the  GB  (CL  or  with  an 
impurity)  is  just  the  one  emulated  by  the  8-atom  supercell  shown  in  Fig.l.  In  order  to  find  the  GB 
structure  corresponding  to  a  minimum  of  energy,  Eq.(2),  the  interplanar  distances  were  varied, 
while  the  interatomic  spacings  and  the  structure  within  the  (111)  planes  were  left  unchanged. 

The  total  energy,  E,  was  calculated  using  a  modified  EAM  approach: 

E  =Ir  Ee,rt,(n(R))  +1/2Zh3'  V(R3’)  +ESE(n(Rimp))  (2) 

where  Ee^b  V(R,R')  are  the  EAM  embedding  energy  and  the  pair  potential  as  found  for  the 

bulk  BCC  Me  (we  used  the  Finnis-Sinclair  functions  and  parameters  for  Fe  and  W  (11)).  The  third 
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term  is  the  energy  of  the  impurity  atom.  R  and  R'  are  the  positions  of  the  host  atoms,  Rimp  is  that  of 
the  impurity,  and  n(R)  and  n(Rin,p)  are  the  electron  charge  densities  at  the  site  of  a  host  atom  and  the 
impurity  respectively.  The  electron  charge  density  at  a  given  site  was  taken  to  be  a  superposition  of 
the  free  atom  charge  densities. 

As  follows  from  our  calailations  the  volume  difference  between  the  CL  GB  and  that  with  an 
impurity  increases  monotonically  with  the  impurity  atomic  number.  From  the  intuitive  point  of 
view,  the  increase  in  volrune  is  expected  to  result  in  a  weakening  of  interatomic  bonding,  though,  in 
principle,  an  impurity  may  exert  a  stronger  interaction  in  spite  of  the  lattice  expansion. 

Fig.  3  schematically  shows  the  atom  arran¬ 
gement  in  the  GB  core.  The  counterparts  of  the 
nearest  neighbor  distances  in  the  bulk  BCC 
lattice  (0.2476nm  in  Fe  and  0.2741nm  in  W)  are 
the  distances  between  nearest  atoms  in  the  [111] 
direction:  Mel  and  Me4  (Fel-Fe4  and  W1-W4  - 
see  Fig.l).  In  the  CL  GBs  these  distances  are 
larger  than  in  the  bulk,  while  the  shortest 
distances  are  the  ones  between  two  Me2  atoms 
>  (Fe2-Fe2  and  W2-W2)  across  the  GBs.  In  the  CL 
GBs  there  are  significant  voids  (occupied  by  an 
empty  sphere,  O  in  our  model);  the  distances 
between  atoms  Me3  across  the  void,  and  thus  the 
GB  plane,  (Fe3-Fe3  and  W3-W3)  are  quite  large: 
0.3108nm  in  Fe  and  0.3614nm  in  W.  Thus,  in  the 
CL  GB  the  strongest  interaction  is  Me2-Me2 
across  the  GB,  followed  by  Mel-Me4  and  Mel-Me2.  With  an  impurity  atom  in  place  of  the  empty 
sphere,  the  interatomic  interaction  fundamentally  changes.  Now  the  shortest  distance  is  that 
between  an  impurity  atom,  •,  and  Me3,  the  distance  between  two  Me2  atoms  across  the  GB  being  the 
second,  l^e  interaction  between  the  impurity  atom  and  atom  Me3  becomes  of  utmost  importance, 
since  this  is  the  interaction  that  is  actually  responsible  for  the  intergranular  cohesion:  with  the  Me2- 
Me2  distances  increasing,  the  main  force  that  keeps  the  two  halves  of  the  crystal  together  is  the  one 
between  the  Me3  atoms  via  the  impurity. 

It  was  also  found  that,  both  in  Fe  (25)  and  W  (26)  GBs,  the  interplanar  distances  oscillate  as  a 
function  of  the  ^stance  from  the  GB,  the  deformation  waves  decaying  by  the  10th- 12th  plane  away 
from  Uie  GB  (plane  1  is  the  GB  plane).  An  interesting  feature  of  the  CL  GB  relaxation  is  that  the 
distance  between  the  2nd  and  3rd  planes  is  only  a  half  of  the  (111)  interplanar  distance  in  bulk  BCC 
Fe  and  W  (0.0407nm  vs  0.0825nm  in  Fe  and  0.0550nm  vs  0.0914nm  in  W).  Though  neither  Fe  nor  W 
undergo  transformation  into  the  ca-phase,  the  "misbalance”  in  interatomic  interactions  arising  due 
to  the  GB  results  in  the  tendency  for  plane  3  to  nearly  collapse  into  plane  2  (the  u-phase 
configuration). 

Both  in  Fe  and  W,  impurity  atoms  result  in  some  "damping”  of  the  relaxation  deformation 
waves,  i.  e.  decreasing  the  oscillation  amplitudes.  This  damping  is  most  pronounced  for  B  and  C. 
Although  the  distance  between  planes  1  and  2  (which  is  half  the  distance  between  two  Me2  atoms 
across  the  GB)  monotonically  increases,  the  tendency  of  plane  3  to  nearly  collapse  into  plane  2 
disappears:  both  in  Fe  and  W,  in  the  progression  B  through  O  the  Me2-Me3  distance  is  almost  equal 
to  that  in  the  bulk.  Except  for  the  Me2-Me3  distance,  the  amplitude  of  the  deformation  wave  increases 
wi^  Al,  Si,  P  and  S.  Damping  the  deformation  wave  may  be  interpreted  as  "cohesion  enhancement", 
while  the  corresponding  increase  of  the  deformation  wave  oscillations  may  be  thought  of  as  resulting 
in  "decohesion”. 

These  conclusions  based  on  semi-empirical  calculations  also  agree  with  our  earlier  first- 
principles  analysis  (8,9,26).  It  was  shown  that  upon  introducing  an  impurity  atom  in  the  GB  the 
impurity's  valence  electrons  hybridize  with  the  Fe  or  W  sp-  and  mostly  d-  electrons.  As  a  result 
covalent  bonds  are  created.  Both  in  Fe  and  W,  B  happens  to  develop  the  strongest  hybridization  with 
the  Me's  electrons.  The  weakening  of  covalent  bonds  •-Me3  (the  most  important  for  intergranular 
cohesion)  with  increase  of  the  impurity's  atomic  number  is  most  pronounced  in  Fe.  H  is  known  to  be 
one  of  tibe  strongest  embrittlers  in  steels.  The  first-principles  calculations  (9)  show  that  while  the  H- 


Fig.  3.  Atom  arrangements  in  GB  cores: 
a)CL  GB;  b)  GB  with  .an  impurity  atom 
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Fel  covalent  bond  (in  the  GB  plane)  is  reasonably  strong,  the  H-Fe3  bond,  across  the  GB,  is  extremely 
weak.  The  hydrogen's  only  electron  does  not  participate  in  metallic  bonding.  Contrary  to  a  gener^ 
belief,  it  is  not  transferred  to  Fe  d-bands,  but  rather  stays  localized  in  the  impurity  band  situated 
below  the  Fe  valence  bands. 


Fig.  4.  AE=Egb(^)*Cgb(CD>  energy  difference  between  the  GB  with  impurity  •  and  CL  GB,  vs 

the  Periodic  Chart  group  number. 

The  analysis  of  GB  energetics  reveals  only  one  aspect  of  an  impurity  embrittling  effect.  From  a 
thermodynamic  point  of  view  (27),  the  impurity's  embrittling  potency  depends  on  the  difference 
between  the  free  energies  of  the  impurity's  segregation  on  the  initial  GB  and  on  the  two  free  surfaces 
emerging  upon  fracture.  The  higher  the  difference,  the  stronger  the  embrittling  potency  of  the 
impurity.  Thus,  in  order  to  have  a  complete  thermodynamic  picture  of  the  embrittlement,  one  should 
be  able  to  analyze  the  enei^etics  of  the  same  impurity  on  the  free  surface.  As  a  much  less  rigorous 
but  simpler  criterion,  in  Ref.(5)  the  sublimation  energy  differences  between  the  host  and  impurity 
atoms  were  calculated  in  an  ideal  solution  model  for  over  60  elements.  According  to  Ref.  (5),  among 
the  elements  analyzed,  only  B,  C  and  Os  may  be  cohesion  enhancers  in  W.  In  Fe,  over  a  dozen 
elements  can  play  that  role.  In  our  more  rigorous  approach,  the  effects  of  impurities  on  GB  stability 
can  be  analyz^  by  simply  comparing  the  GB  energy  differences,  AE,  between  the  GB  with  impurities 
and  the  CL  GB. 

The  corresponding  values  for  the  impurities  discussed  are  plotted  in  Fig.  4.  One  can  see  that 
both  in  Fe  and  W,  the  GB  stability  decreases  from  B  towards  O,  and  the  energy  difference  becomes 
positive  for  Si,  P  and  S  (in  Fe),  and  P  and  S  (in  W)  --the  strongest  embrittlers.  The  latter  means  that 
GBs  with  the  corresponding  impurities  make  the  (IBs  unstable  at  absolute  zero. 

CLaaclttgion 

The  site-competition  effect  is  well  known  and  experimentaUy  observed  in  many  systems.  In 
this  paper  we  made  an  attempt  to  understand  this  phenomenon  on  the  microscopic,  electron-atom 
level.  From  the  theoretical  point  of  view,  the  site-competition  effect  can  be  used  in  order  to  cleanse  the 
GBs  of  the  "harmful"  impurities,  such  as  P  and  S.  Boron  may  be  such  a  cleansing  agent.  Not  only 
does  it  tend  to  displace  other  metalloid  impurities  from  the  GBs,  it  was  also  fbimd  to  improve 
intergranular  cohesion.  This  effect  of  boron  has  been  observed  in  other  systems  also  (26).  An 
advantage  of  a  "site  competition"  cleansing,  as  compared  with  traditional  gettering,  would  be  the 
absence  of  a  second  phase  precipitate,  which  may  be  difficult  to  control.  On  the  other  hand,  to 
"deliver"  B  to  the  GBs  to  be  dean^  may  be  quite  a  difficult  problem.  In  our  analysis,  we  completely 
disregarded  alternative  chemical  reactions  of  B  both  with  the  impurity  and  the  matrix  atoms  that 
might  result  in  formation  of  boron  compounds,  such  as  metal  borides,  boron  oxides,  etc.  The 
thermodynamic  analysis  of  the  possibilities  of  such  alternative  reactions  would  be  very  helpful.  We 
believe  that  an  experimental  investigation  of  the  effect  of  B  on  ductility  of  Fe  and  W  is  of  extreme 
importance.  Recently  (15,16),  it  was  shown  that  B  introduced  in  ppm  quantities  into  an  Fe-P  alloy 
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significantly  improves  the  alloy  ductility.  At  the  same  time,  the  Auger  analysis  showed  the  decrease 
in  the  P  concentration  on  the  fracture  surfaces  -  a  manifestation  of  the  site  competition  effect. 
Provided  the  technological  difiicxilties  can  be  overcome,  microalloying  with  B  could  result  in 
essential  toughness  improvements  in  current  existing  high-strength  steels.  10-50  ppm  of  B  can  also 
make  W  less  brittle.  Experimental  research  on  microalloying  W  with  B  is  currently  in  progress  at 
Army  Research  Laboratory  in  Watertown,  MA. 
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